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PRESENT STATUS OF LITHIUM-DIFFUSED 


SILICON SOLAR CELLS 


I. INTRODUCTION 

The present report is concerned with lithium-diffused silicon solar cells. 

First, a brief history of work on lithium-diffused semiconducting silicon will be 
given, and a few experimental results on the radiation effects and the after- 
effects will be discussed. Second, the radiation damage and recovery of 
lithium -diffused solar cells will be analyzed from the physical point of view. 
Finally, experimental results of the initial efficiency, and the subsequent radi- 
ation damage recovery will be discussed. 

The solar cells used in this study were made by RCA Laboratories, Heliotek, 
Hoffman Electronics Corp. , and Texas Instruments, Inc. The enthusiasm and 
active cooperation of these suppliers are appreciated. The experiments were 
performed by Messrs. G. Meszaros, Y. M. Liu, and W. McGunigal. The dis- 
cussions among the members of the Radiation Physics Section of Goddard Space 
Flight Center were of great help in planning and carrying out the program. In 
addition, encouragement of Mr. M. Schach on this project is hereby acknow- 
ledged. 


1 



n. HISTORICAL REVIEW 

Diffusion of lithium (Li) in silicon (Si) was first investigated by Fuller and 
Ditzenberger (Reference 1) by the p-n junction method, in which the location 
of the junction, which is formed by diffusing Li into p type Si, is measured. 

They found that 

D = 94 x 10 -3 exp ^-9. Ill cm 8 sec -1 , (II - 1) 

D is the diffusion coefficient and T is the absolute temperature. This leads to 
a room-temperature value of D = 1. 3 x 10 -14 cm 8 sec -1 . At 400° C, D will be 
3.58 x 10“ 9 cm 8 sec -1 . 

In a later work (Reference 2), Fuller and Severiens gave, by measuring mobility 
vs applied field, 

D = 2. 3 x 10 -3 exp^-7. 65 x cm 8 sec -1 . (II - 2) 

This would give D = 2. 0 x 10“ 14 cm 8 sec -1 at room temperature and 
2.5 x 10 -8 cm a sec -1 at 400° C. In this experiment, Li is again diffused into 
p-type Si. Because there is an interaction and precipitation of Li in p-type, 
Boron (B) doped Si (see Reference 3), the value of D might be different from 
that obtained when Li is diffused into n-type Si or pure Si, but no data pertaining 
to these cases are available. Furthermore, since Li interacts effectively with 
oxygen (O) to form LiO + or other forms of the oxide, which are stable at room 
temperature (Reference 4), it should be noted that all these diffusion studies 
were carried out with Si crystals grown from a quartz crucible, i. e. , Si of high 
O content. A recent work of Gilmer, et al. (Reference 5) gives an up-to-date 
and complete study of Li-O complexes as donor impurities in Si. 
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The investigation of the existence of the interaction of Li with the complex 
defects introduced by irradication was initiated by Vavilov and his co-workers 
(Summarized in Reference 6). By observing the change of the characteristics 
of the temperature dependence of carrier concentration, and of the lifetime of 
minority carriers in p-type Si, with or without Li doping, the following obser- 
vations were made: 

1. hi n-type Si which was obtained by inversion from p-type Si, radiation 
damage results in reversion to the original p-type. This is explained as due to 
the interaction of Li with the radiation-induced defect. The donor concentration 
contributed by Li is reduced. This effect is called precipitation by Vavilov. 

2. The rate of production of A -centers (O vacancy complex) in Si contain- 
ing Li is less than that in Si containing phosphorous (P) or arsenic (As). This 
is attributed to the formation of the LiO + complex. Thus, the effective oxygen 
concentration which provides a component of the A-center is reduced. 

3. The degradation of the lifetime due to electron irradiation in Si con- 
taining Li shows a saturation for 1 Mev electrons after about 1 X 10 15 e/cm 3 . 

In Si containing P but no Li, no such saturation, indeed, only a continuous 
degradation, is observed. 

We will discuss more specifically the first two aspects later. The third aspect 
has not been verified (Reference 7). Although it was not expressly stated, pre- 
sumably all crystals studied by Valvilov et al. have a high oxygen content. 

Vavilov attributes the idea that the disapperance of defect centers is due to the 
interaction of impurity centers with structural defects, to Roberts et al. (Refer- 
ence 8). In their brief paper, Roberts et al. found an increase in the majority 
carrier lifetime with an increase of the oxygen concentration, in contradiction 
to the usual expectation that any departure from crystal perfection will reduce 
the lifetime. These authors presented an interesting hypothesis that there is an 
unknown center which either can be nullified by chemical interaction with oxygen 
or competes with oxygen for a place in the lattice. This point has not been 
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further verified. In the meantime, a successful application of this idea to the 
case of defects introduced by irradiation is clearly shown by Vavilov et al. 

It was the lecture of Vavilov (Reference 6) that revived the interest of 
Rappaport in investigating the introduction of Li into Si solar cells. * Immedi- 
ately a study was made of the surface barrier diode with Li diffusion (Reference 
9). The degradation of the diffusion length was measured as a function of the 
radiation flux. A recovery of the diffusion length was observed after cessation 
of the radiation and a waiting period of several hours. 

An attempt was made to apply this idea to solar cells with the following configur- 
ations: (1) n-on-p solar cells, (2) p-on-n solar cells made from pulled crystals, 
and (3) p-on-n solar cells made from float zone crystals but with a high P con- 
tent (~1 Q cm or less). In all cases, results the same as those observed in 
the case of the surface barrier diode were not obtained. 

It was the contribution of Wysocki, at this point, to experiment on p-on-n solar 
cells made from float zone crystals with a low P content (larger than 10 cm). 

The result was spectacular: For the first time, a practically complete recovery 
at room temperature was observed in the solar-cell performance which had been 
degraded by the radiation. For 1 Mev electron irradiation and for an accumu- 
lated flux of 1 x 10 14 e/cm 3 , a complete recovery was observed in a matter of 
hours. If the flux is higher than this value, recovery is less complete. This 
phenomenon is also dependent on the rate of radiation flux. The above irradia- 
tion was carried out at a rate of 1 x 10 14 to 1 x 10 16 e/cm a hr. If the flux rate 
is 5 X 10 13 e/cm 8 hr in fact, because the rate of the recovery is compatible with 
the rate of damage introduction, no damage is observed at all. All these obser- 
vations have an analogy in the annealing of n-on-p solar cells (Reference 10). 

The critical flux levels are equal in both cases, within the accuracy of the exper- 
iment. If recovery is permitted intermediately by reducing the radiation rate, 


♦Private communication. 
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the recoveries are complete in both cases. In fact, in this way the total flux 
can exceed the critical flux obtained the result of a fast irradiation rate. 

This brings the status up to early 1966. The progress made since is the main 
material of the present report. 


EL PHYSICAL MODELS OF SPONTANEOUS ANNEALING 
The present models of the mechanism of the spontaneous annealing of Li- 
diffused solar cells are based on indirect observations and are mostly specula- 
tive. Two different models exist, one in which the role of Li is to compensate 
the defects (Reference 11) and one in which the role of Li is to form metastable 
defects (Reference 12). 

The Role of Li is to Compensate the Defects 

When the Li is to compensate the defects, Li + is moving randomly in the 
crystal lattices by diffusion. When a defect due to the radiation occurs near 
the instantaneous site of Li + , a defect representing a recombination or trapping 
state to the minority carriers in the n-region would have a negative charge; 
therefore, an interaction with Li + would lead to a neutralization of the charge 
state and a consequent annealing. 

One support of this model is the experimental data of the temperature depen- 
dence of the recovery time (Fig. El-1). The isothermal recovery times r at 0°, 
25°, and 78° C are measured at toom temperature. For this limited tempera- 
ture range and radiation fluxes of 1 and 3 x 10 14 e/cm a of 1 Mev electron radia- 
tion the plot of t vs. l/T, where T is the absolute temperature, is a straight 
line whose slope of 0. 66 ev corresponds closely to the value of the activation 
energy of diffusion previously reported (Reference 1 and 2). But in view of the 
lack of the precise value of the activation energy of diffusion and the limited 
temperature data used, this cannot be treated as conclusive evidence. There 
is also some evidence that, even if the annealing is controlled by the diffusion 
of Li to the compensation site, a complexity could result from the recent obser- 
vation that the recovery rate is dependent on the density of the defect. * This 
is not taken into account in the treatment of the diffusion theory from which the 
experimental results were analyzed. 

*P. H. Fang and Y. M. Liu (to be published) . 
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Probably the strongest argument against this model is that the center which is 
assumed to be annealed by Li cannot be identified with any presently known 
defect. * As described in the beginning, the spontaneous annealing is observed 
in float zone grown, but not pulled, Si; this excludes the importance of the A- 
center. Secondly, the annealing is observed in moderatc-to-low P concentra- 
tion material in, and low concentration B material which undergoes subsequent 
conversion to n-type Si; thus the recovery is independent of E-center defects. 
Finally, the divacancy defect cannot be important because of its low rate of 
production. 

Another difficulty of this model is the following observation; Solar cells, after 
electron irradiation and after an almost complete recovery, are heat-treated up 
to 350° C for several hours to look for signs of reverse annealing. In this tem- 
perature range, all known defects are annealed by a dissociation of the vacancy 
from the complex. * In the case of Li, therefore, a dissociation of Li from the 
defect, or a dissociation of the defect itself, should lead to an observation of a 
change in photovoltaic characteristics. The result is negative — no ascertain- 
able changes were observed. 

This deductive reasoning suggests the simpler and more direct model which 
will be discussed next. 

The Role of Li is to Form Metastable Defects 

When high-energy particles enter the Si lattice, vacancy-interstital pairs 
will be created. Both vacancies and interstitials are highly unstable at room 
temperature, but vacancies are less unstable and most of them are transformed 
into stable states by forming complexes with impurities or with other vacancies. 
In float zone Si with Li as the major impurity, there is a great probability that 
there will be an interaction between the vacancy and Li. This interaction has 
been observed by Vavilov (Reference (6). Goldstein has conjectured the exist- 
ance of a Li-V (where V denotes vacancy) complex through electron paramag- 
netic resonance studies (Reference 13). 

*G. Watkins, loc. c. f. Referenced, p. 97. 
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The model of the Li-V complex is an interstitial Li + and a V” in a lattice posi- 
tion (Reference 12). Since Li + forms a stable state at an interstitial position 
in the regular Si lattice, even under Coulomb attraction by one of the four 
nearest neighbor Si ions, the Li + probably could maintain a metastable state. 
The annealing mechanism in this model, the activation energy of the recovery 
time in the microscopic diffusion, is assured to be approximately equal to that 
of the macroscopic diffusion. 

Much basic information is needed to determine the more valid of the present 
phenomenological models. This knowledge is equally needed in the evaluation 
and the projection of the present result to technical applications. In this re- 
spect, an important question is that of stability over extensive periods of time 
and over a wide range of temperature. Also a very interesting question is 
whether the recoverability from a large desage of radiation flux is limited by 
the initial Li concentration. If each annealing requires one Li + ion in the crys- 
tal, eventually, when the number of vacancies becomes approximately equal to 
the number of Li + ions, the spontaneous recovery will have ceased. As will be 
discussed later, we have noted that the defect production rate as well as the 
degree of recovery depends appreciably on the concentration of the instantan- 
eous defect density, and, under some dynamic conditions, it seems that the 
extent of recovery after an extensive accumulated flux dose does not seem to 
be limited by the available Li concentration. A possible mechanism will now 
be given. 

According to the second model of annealing, the recovery of the radiation dam- 
age is due to the occupancy of the lattice position by Li + of the Li + -V“ complex 
defect. This position requires the trapping of an electron by Li + , thus, 

Li + + e -* Li° , 
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is thermodynamically not stable. The most stable configuration is 

Li^ + Si. -♦ Li.+Si. 

1 1 "b 

+ 

-* Li. + Si. + e, 

1 </ 

which describes an interchange between interstitial (i) and lattice (l) positions 
and a subsequent capture of a hole. Two sources of Si. are 

1. A Frankel defect of thermal origin, and 

2. In the case of simultaneous irradiation, continuously created Si. 


The existence of such processes will be vital in determining the ultimate 
limitation of the practical use of Li solar cells. 


10 



IV. DIFFERENT LI SOLAR CELLS 

On the basis of the present knowledge of the interesting properties of the inter- 
action of the radiation-induced defect with Li, four types of Li solar cells have 
been developed: 


Type I. 
Type n. 
Type m. 
Type IV. 


p-on-n with float zone Si, 

p-on-n (phosphorus) on n (lithium) with float zone Si, 
p-on-n with pulled Si, and 

n-on-p with float zone Si, with a relatively small Li concen- 
tration in the p -region. 


Type I Cell . The Type I cell is shown in Figure IV-1. This is the prin- 
cipal type of Li cell. The important constituent of the cell is a Si crystal of 
low O content, a p (up to now, boron) -type front region, and Li as the domi- 
nant donor in the n-region. Most discussions up to now were based on this 
type of solar cell, and it will be described exclusively in later sections. Since 
a sufficiently specific description of the cell without misnomer requires a cum- 
bersome expression and since the cell will be frequently referred to in the 
future, it seems most desirable to coin a special name for this type of cell. 

It is remarkable that the cells were made by Wysocki before we knew or had 
any clear physical picture of the nature of the radiation damage. The foresight 
of Wysocki (Reference 7) was instrumental in the cells discovery and it seems 
only proper to designate this type of cell as the Wysocki cell. Therefore, the 
type I cell will be so designated hereafter. 


Type II Cell. The principle of the type II cell is based on the following 
observation: Because of a compensation near the p-n junction and because of 
the diffusion profile configuration, there is a depletion of Li concentration in 
the n-region near the p-n junction. With irradiation, further reduction of free 
Li results, because of a complex formation between Li and the radiation-induced 


11 



2 CM- 



12 


I -A 


defect. This diminishing of Li will be most effective near the junction because 
of the original comparatively lower Li content in this region. 

An approach to solving this problem is to make double n- layers*. The first 
n-layer, consisting of the p-n junction and the immediate n-region, is a few 
microns thick and is made by the diffusion of P into Si, and the layer thickness 

is controlled either by a diffusion or by a subsequent lapping. The p-layer is 
then diffused into the phosphorus n-layer to form the p-n junction. A Li dif- 
fusion, in the same manner as in the Wysocki cell, is then performed. 

Figure IV-2 is a diagram of a Li Type n cell. The radiation damage property 
of the type n solar cells will be discussed in the Section entitled "Some Experi- 
mental Results Pertaining to Space Applications." 

Type in Cell When Si of high O content, as is found in pulled crystal, is 
used, the defect induced by the radiation is stable at room temperature. How- 
ever, an effective annealing stage occurs at a temperature of 120° to 150° C 
(Figure IV-3.) This is to be compared with the effective annealing tempera- 
ture of about 400° C in the usual n-on-p solar cell, t The nature of this defect 
has been conjectured as a vacancy forming a complex with the Li-0 pair. 

From a practical point of view, if a high- temperature annealing is to be used 
as a means of solving the problem of the radiation damage, a 150°C annealing 
is much to be preferred from the point of view of both the demand of the stabil- 
ity of the components and connections in the annealing machinery and in the 
required thermal power of annealing. 


*P. H. Fang, NASA Goddard Space Flight Center (Patent pending) . 
tP. H. Fang, Astron. Acta (to appear). 
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Type IV Cell. In the original work of Vavilov (Reference (!) an interesting 
observation was reported: In p-type (B-doped) Si (Reference when Li is 
introduced, an over-compensation resulted in the conversion to n- type Si. A 
subsequent electron irradiation showed a reversion to p-type. The interpreta- 
tion was that Li-B pairs are formed and the effective original acceptors are 
reduced in number. With electron irradiation, the defects formed (presumably 
single vacancies) interact strongly with Li, B is reliberated, and the effect 
results in a restoration of the majority carrier concentration. 

Suppose now a n-on-p solar cell is made from p-type Si with a high concentra- 
tion of B, say, corresponding to 1 fi-cm. If the introduction of Li resulted in 
the formation of p-type Si of 10 fl-cm, the solar cell is now a 10 O-cm n-on-p 
cell. After irradiation of this type of cell, assuming there is no radiation effect 
other than the reconversion of high resistivity p-type to low resistivity p-type, 
a radiation-induced improvement in the solar-cell performance could be 
expected. 
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V. LI SOLAR CELL TECHNOLOGY 

The preparation of Li solar cells has been described previously in the work of 
Topfer et al (Reference 14). More up-to-date reports on the technological pro- 
cesses will become available in the near future (Reference 15). Therefore, the 
description here will emphasize the physical aspects and will leave out the tech- 
nical processing. 

Figure IV-1 shows the schematic of a p/n, Li-diffused, Si solar cell. In the 
fabrication, the only difference from the ordinary p/n solar cell of B (p-type) 
and P (n-type), is that a Li diffusion is carried out from the base of the solar 
cell. Since the diffusion coefficient of Li is enormously greater than that of the 
presently used p-layer dopant (B), the Li diffusion is carried out after the high 
temperature process of p-n junction formation by B diffusion. The Li diffusion 
is then carried out near 400° C. 

There are three commonly used methods of diffusing Li into Si: 

1. Paint the Si with a Li-oil emulsion. This is the simplest method, but 
the amount of Li applied to the surface cannot be easily controlled. The diffu- 
sion is carried out in a high-temperature furnace. 

2. Immerse the Si in a Li-Sn alloy bath. The temperature of the alloy in 
the liquid state is about 400° C. The diffusion time required is usually longer 
than in 1 , but a better uniformity can result. This method also presents some 
difficulty in the protection of the p-layer and preventing edge leakage. 

3. Evaporate Li into Si with a subsequent diffusion. This method allows 
for better control of the amount of Li applied to the Si surface. In comparison 
with the Li-oil emulsion, since the oil is much less pure, the Li evaporation 
directly from Li metal reduces the introduction of unknown or undesired 
impurities. 
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On the problem of the solubility of Li in Si, the following parameters are prob- 
ably important and are presently largely unknown. 

1. The optimum concentration of Li. From the radiation damage and 
annealing point of view, high Li content is desirable. But high Li content could 
also result in a larger carrier removal rate, as observed in an analogous case 
in p-type Si (Reference 16). 

According to the experimental data of Pell (Reference 17), the solubility of Li 
at 27° C is 10 13 cm -3 . The concentration of Li obtained by the current diffusion 

X 6 1 7 — 3 

practice is of the order of 10 to 10 cm . Therefore, the system is highly 
super-saturated. The consequence is an expectation of the precipitation. The 
quantitative value of the precipitation time is unknown at the present. 

2. The effect of the impurity on the solubility of Li. There is an obvious 
difference in the case of p-type Si from that of n-type Si. This is qualitatively 
explained by the electrolytic solution model (Reference 18). The Si specimens 
used in earlier experiments were pulled crystal of impurity concentration, 
before the Li diffusion, of the order of 10 14 to 10 le impurities cm -3 . Since in 
the pulled crystal the O content is greater than 10 16 cm” 3 and there is a strong 
interaction between Li and O to form LiO” and Li a O, the problem of what will 
be the solubility of Li in float zone Si where the O content will be two or more 
orders of magnitized less than that in pulled crystal is still to be investigated. 

Another problem is the solubility of Li in very pure Si of impurity content of 
the order of 10 l3 cm -3 , corresponding to a resistivity of about 10 3 fi-cm. If 
one is to emphasize the room temperature annealing, any impurity, other than 
Li, should be minimized because, according to present knowledge, the defects 
associated with all impurities, except those of Li, are stable up to 150° C or 
higher. Therefore, it is natural to start with Si material with as low an im- 
purity concentration as possible. 
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3. The effect of dislocations on the solubility of Li. Two kinds of crys- 
tals used currently in Li solar cells are float zone and Lopex. * While float 
zone material has a dislocation concentration in the range of 10 14 cm -3 , Lopex 
has about two orders of magnitude less. There is some indication that the dif- 
fusion coeficient of Li in Lopex material is smaller than in float zone 
material. ** 

4. In connection with the solubility, there is also an interesting problem 
of the concentration profile. Since in the Li solar cell the donor is introduced 
by diffusion, and, furthermore, since a p-n junction is present during this 
diffusion, the region having an excessively high concentration, greater then 
10 19 cm“ 3 , of acceptors behaves as an infinite sink. Therefore, the concen- 
tration of Li cannot approach a constant value throughout the depth of the cell. 
This unintentional concentration profile actually satisfies the configuration of 

a p-on-n drift field solar cell and was considered by Ross long before its radi- 
ation and annealing properties were investigated. *** 


*A product of Texas Instruments, file. 

**D. Kendall (private communication). 

***B. Ross (assigned to Hoffman Electronics Corp.). Patent pending. 
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VL RADIATION DAMAGE AND RECOVERY 

The general setup for studying the photovoltaic performance of solar cells before 
and after radiation, and the history of recovery is similar in practice to that of 
the early studies of the usual n-on-p solar cells except in one aspect: In a n-on- 
p case, only the initial characteristics and the radiation damage are measured; 
the recovery, if any, will be much less than 1 percent. 

In our measurement, less emphasis is placed on the assimilation of space con- 
ditions both with respect to radiation and to photo-illumination, and more 
emphasis is placed on a simple and repeatable measurement. In view of the 
great gap existing between the basic physical properties and their relation to 
device parameters, our philosophy has been to look at only the gross device 
characteristics. Specific basic studies are made only when some significant 
effects are recognized. In this way, information for practical purposes is pro- 
vided. At the same time, in recognition of the great complexity of the elec- 
tronic device — the solar cell in the present case — as a physical system, 
unless the effect of degradation of the recovery is large, it is almost meaning- 
less to attempt to describe an effect as significant, which might be extraneous. 

The equipment and experimental routine for the above investigation are quite 
simple and have been described elsewhere (Reference 19). The task in all the 
measurements is to compare solar cells made by different processes or irradi- 
ated under different conditions. 

The primary investigation of the radiation damage in the Wysocki cell uses 
1 Mev electron irradiation at room temperature. Unless otherwise specified, 
this will always be the condition of radiation. When the current-voltage (i-v) 
curve is measured for the cell before and after irradiation, typical data like 
those of Figure VI-1 are obtained. The i-v curves in the figure represent 
isothermal recovery from an initial damage which is shown as dotted lines: 
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For the practical reason that the location of the radiation source and the site 
where the measurements were made is not contiguous, the points in the dotted 
line, were obtained by extrapolating the isothermal data back to zero time. 

From this, we estimated the initial damage as of the order of 20-pereent after 
a 1 Mev, 10 14 e/cnr ! irradiation. This damage is comparable to that of the 
ordinary p-on-n (no Li) solar cell and is considerably heavier than that of the 
ordinai’y n-on-p solar cell after a similar irradiation. 

The characteristic time in the recovery process is on the order of hours. 

Since the rate of electron irradiation is 1 x 10 le e/cm 2 hr, the radiation time 
for an accumulated flux of 1 x 10 14 e/cm s , which is less than 1 minute, is 
negligibly in comparison with the recovery time. When we study the case of 
large flux damage, more complications arise, and additional considerations 
will be made. 

Therefore, in the first study we will confine ourselves to low flux levels of 
electron radiation. In the development of the analysis, more specific aspects 
such as the dependence on the temperature, the flux level, the flux rate, and 
the defect state will be discussed separately. 

Solar cells made by the best methods presently available, with some variation 
of the temperature and the time of diffusion were received from the contractors 
(Reference 15) at different times and irradiated with 3 xl0 14 e/cm s of 1 Mev 
electrons. Recoveries were measured periodically up to several weeks in some 
cases. From the data, we obtained histograms (Figures VI-2, VI-3, and VI-4) 
of the relation between the initial performance characteristics and the degrad- 
ation due to the radiation of 115 solar cells. This statistical result is prelim- 
inary. At a later stage, when processes can be better controlled, a more 
specific analysis will be made. The histogram shows the value after degrada- 
tion and recovery as a function of the original value of the short circuit current 
(i), open circuit voltage (v), and efficiency (77). 
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There is a slight correlation of a general increase of unrecoverable damage as 
the photovoltaic parameters of i, v and n are increased. But there is also evi- 
dence of a considerable scattering of data which leaves a question unanswered: 
Whether the correlation is an intrinsic property of the Wysocki solar cell, or 
whether it merely represents a transitory technological problem. This ques- 
tion may soon be answered when large numbers of cells are produced under 
better controlled conditions. At present, the physical evidence does not pro- 
vide a definitive correlation, at least in the low flux level case considered here. 
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VH. DEFECT CONCENTRATION DEPENDENCE OF DEFECT 
INTRODUCTION RATE 

The study of the defect introduction rate in the Wysocki cell is complicated by 
the simultaneous occurrence of annealing. Two simplified experiments can be 
made: 


1. The rate of defect introduction is much greater than the rate of 
defect annealing, and 

2. The reverse of 1. 

The van de Graaff irradiation, at a rate of 10 16 to 10 16 e/cm 3 -hr will approxi- 
mate the first experiment when a total flux of less than 1 x 10 lS e/cm® is used. 
Beyond this total flux level, the introduction and the annealing rate become 
comparable. On the other hand, in ordinary space radiation conditions the 
situation is closer to that of the second experiment. While the actual space 
radiation experiment is to be carried out in the near future, another experiment 
which will be closer to the space radiation one is also in preparation, which 
will provide a prolonged period of continuous radiation by means of a y-source 
or isotope source. This investigation is in progress. 

Figure VII-1 shows the radiation damage of 10 cells with initial spread also 
indicated. Sometimes the practice in the literature has been to use a single 
cell and irradiate it to successive high flux levels. Because there is evidence 
that the damage rate depends on the flux level as well as on the damaged state 
of the cell, fresh cells are used in each flux level experiment. The data are 
analyzed with the usual Equation of the degradation of the lifetime, 

- = - + K§, VB-1 

T T o 
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where t and r are the lifetime before the irradiation and after the irradiation 
o 

and recovery, respectively, $ is the flux in e/cm 3 , and K is an adjustable 
parameter. In this equation, a shift of Fermi level due to the carrier removal 
is neglected. Since the diffusion coefficient D and diffusion length L are 
related by 


D = L 3 , vn-2 

T 

and L is proportional to the short circuit current i, we have 

-|=“ VH-3 
i i 

o 

By fitting data at $ = 1 x 10 15 , 

K' = 1. 51 X 10" 13 (cm 4 /amp 3 -e). VII-4 

Equation VII-3 is computed with this K 1 value and is shown as a continuous 
curve in Figure VII-1. The curve fitting is excellent in the low flux region. 
The deviation at 3 x 10 15 e/cm 3 can be explained by the spontaneous recovery 
taking place during the irradiation, that is, the initial damage after the 
cessation of the radiation is lower than that which would be caused by a pulse 
of 1 x 10 15 e/cm 3 . 

The same observation applies to the data of 1 x 10 16 e/cm 2 radiation. In 
this case, an extrapolation from the data of 3 x 10 15 would indicate a stronger 
deviation from the value predicted by Equation VII-3. The weaker deviation 
observed can be explained by an onset of a steeper damage introduction rate 
near this high flux. 
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VIE. DEFECT CONCENTRATION DEPENDENCE OF ANNEALING RATE 
In the study of the annealing of radiation- induced defects in semiconductors 
some evidence was obtained on the dependence of the annealing rate on the 
concentration of impurities (Reference 20). 

When the irradiation is carried out at room temperature, most of the defects 
formed in Si are stable up to several hundred degrees centrigrade*. Therefore, 
isothermal annealing, on which the annealing rate is measured, has to be 
carried out at high temperatures. This introduces some complexities in the 
interpretation of the data since the activation energy of the annealing itself 
could depend on the temperature (Reference 21). This difficulty is minimized 
in the case of annealing of the defects in the Wysocki cell because in this case 
both the radiation and annealing occur at room temperature. Therefore, the 
entire experiment can be carried out at a single temperature. 

Solar cells are irradiated with different total fluxes of 1 Mev electrons at a 
rate of 3 x 10 ia e/cm 2 sec. The specimen is placed on a water-cooled 
aluminum plate, and the temperature rise during the irradiation is less than 
5° C. The isothermal change of the short circuit current of the solar cell at 
room temperature is measured. The data are shown in Figure VIII- 1. 

In the initial stage of recovery, an exponential time dependence adequately 
describes the experiment result. We have, for the fraction of residual 
damage f, 

i - i % exp (-t /t), vni-1 

o a 


* R. R. Hasiguti and S. Ishino, loc. cit. Reference 6; p. 259. 
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Figure VIII- 1 



Where i is the short circuit current, the subscripts o and a stand for initial 
and after annealing, respectively, and t is the time variable. The results 


are given in the following table, 
linear portion of the curves. 

The recovery time is evaluated from the 
Table VEI-l 

Specimen i 0 

ma/2 cm 3 

Total flux/ cm 3 

T x 10 4 sec ma/2 cm 3 

1 

53. 5 

1 x 10 14 

1.1 

2 

53.2 

3 x 10 14 

1.4 

3 

49.4 

3 x 10 14 

1.9 

4 

49.4 

1 X 10 16 

2.5 

5 

50.8 

1 x 10 15 

2.8 

6 

50. 0 

3 x 10 15 

10.6 

7 

49.2 

3 X 10 15 

10.6 


In a more extensive analysis, it is necessary to minimize the contribution 
of p-region degradation. This is achieved by injecting a penetrating light of 
8420 X (with 100.& bandwidth). Furthermore, a linear proportionality between 
the diffusion length and the short circuit current, is assumed and the life time 
of the minority carriers to be proportioned to L -a therefore, the fraction of 
the annealed defects is defined as 

f = (i " a - i " a ) / (i. - i " a ), Vin-2 

'a o f o 

where i, is the short circuit current after irradiation with a total flux f . 
f 

In the computation of f, we have subtracted a portion of the photovoltaic current 

i from i . The current i represents the defference between i and i mea- 
c o c o a 

sured after a very long time when no further change in i is observed. 

3 . 
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This portion of i obviously is contributed by other types of defects and 
therefore is subtracted in our analysis. 

Figure VIII-2 shows the isothermal annealing data of f for specimens with 
different total fluxes of irradiation. We observe a shift of the isothermal 
annealing curve towards the longer time. This corresponds to an increase 
of the recovery time as the defect concentration is increased. 

The important implication of this result is that the dependence of the recovery 
rate on the defect concentration is opposite to the prediction of current models 
which can be collectively represented by (Reference 22) 

f - i+A,T ■ vm - 3 

where is the concentration of defects and X is a positive parameter inde- 
pendent of . From the dimensional structure, therefore (XN^) _1 is a measure 
of the recovery time. Since is a positive increasing function of i. Equation 
VIII-2 predicts a decrease of the recovery time with respect to the concentra- 
tion of defects. 

Taking an arbitrary level of f, we obtain the concentration dependent recovery 
time T which is given in Figure VIII-3. For values of r at f = 0. 50, a good 
approximation is obtained by 

r A _ = 74. 5 exp (2. 86 X 10" 5 i V 3 ). VDI-4 

u. 5 

On the other hand, for values of rat f = 0. 10, 

T o i = 5 ' 06 x 103 exp 73 x 10 ~ 6 VHI-5 
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It is interesting to note that the main difference between Equations VIII-4 and 
Vni-5 is the pre-exponential factor. The numerical values given in Equations 
Vm-4 and VIII- 5 could conceivably be dependent upon the processes and the 
detailed technology of the solar cell fabrication. The importance is that a 
§ 1 / 3 dependence describes well the general behavior. 

From the same data, we can analyze the dependence of the fraction of the 
residual damage on the concentration of the defect. Figure VIII-4 shows the 
result. The curve is calculated from the following equation: 

f= 6.38 X 10" 3 exp (2.12 X 10" 5 I 1 / 3 ). VIII- 6 

In the same figure, 6 = i~ 2 - i “ 2 is also shown. The same f 1 ^ 3 dependence 

$ o 

is also found to be a satisfactory description of the experimental data. 

Within the experimental accuracy, the coefficients in the exponential can 
be treated as if they were about the same. 

We will return to this discussion in the section entitled "Kinetics of 
Annealing. " 
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IX. KINETICS OF ANNEALING 

In the early investigation of the annealing, the experimental data used in making 
comparisons with physical models was exclusively that obtained from isothermal 
and isochronal annealing. (Reference 23). Our investigation shows that the flux 
dependence of the annealing provides a more stingent test of the physical model. 
By this test, we have proved that currently existing annealing models predict 
annealing properties in contradiction to the experimental observations. A 
phenomenological model is proposed* in which an interaction is assumed be- 
tween defects, and consequently the activation energy of defect formation and 
annihilation becomes a function of the defect concentration. 

Since the above work has already been reported in the literature, a repeated 
description, even if it could have some pedagogical value in constructing some 
mathematical approved — the original papers are necessarily brief because they 
appear as letters — will not be attempted here. At this stage, a more impor- 
tant problem is to look for the physical foundation of this model. 

First, we will estimate the activation energy represented by Equation VIII-5: 

rO.l = 5.06 X 10 a exp 2.73 X 10" 6 * 1 / 3 S ec. EX-1 

Since the experiment pertains to room temperature (300° K), we have 

2.73 X 10“ B $ V 3 = E x , IX- 2 

where E 1 is the defect concentration dependent part of the activation energy. 


*P. H. Fang, Phys. Letters . 24 A, (251) 1967. 


40 



In this sense, 


T q - 5. 06 X 10 2 exp (-E 0 /kT) , IX-3 

where E 0 is independent of concentration. From the discussion in the section 
entitled "Li-Solar Cell Technology," we will tentatively assign 

E 0 = 0. 66 ev. IX-4 

Therefore, r 0 = 1 . 01 x 10 -9 sec . On the other hand, 

Ei = 7.06 X 10 -7 $ V 3 IX- 5 

The defect production rate in n-type Si is on the order of 0. 1 to 0. 01. In this 
case, 

Ei = (1. 52 to 3.29) X 10 -7 n 1 / 3 , IX-6 

where n is the defect density. Therefore, for n = 10 1S , Ei = 1. 52 to 3. 29 mv; 
for n = 10 15 , E = 15. 2 to 32.9 mv. This is quite small compared with the E 0 
value of 660 mv; therefore, the E-value will dominate in the measurement of 
the activation energy. 

In the search for a possible mechanism, we note that in the case of Hall con- 
ductivity in germanium (Ge) an impurity concentration dependent activation 
energy is given by (Reference 24, 25), 

E = E 0 - (2. 35 to 4. 3) X 10 _8 n 1 / 3 ev. IX-7 

The numbers involved are not very different. The positive sign in the case of 
Equation IX-6, versus the negative sign in the case of Equation IX-7, for the 
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term n 1 / ; bould be attributed to the increase or the relief of the strains in the 
lattice respectively. 


Finally, according to the previous work*, we have the following kinetic equa- 
tion: 

H? = W (-Ei/kT). IX-8 

After integration and taking n = no at t = 0, 

(a n 0 V 3 ) (o' n 1//s ) = t/1506. IX-9 

This will be the isothermal equation. The value of aia given by Equation IX-6 
and is dependent on the temperature. 


*P. H. Fang, Phys. Letters 24 A, (251) 1967. 
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X. DYNAMICS OF ANNEALING 


The problem to be discussed is the situation one would encounter in space 
applications. In view of the scarcity and the tentativeness of the available data, 
the point to be stressed is the qualitative characteristics rather than the actual 
numerical value of the computations. 

Based on the equation of the kinetics of annealing, if there is a simultaneous 
recovery in the presence of a radiation source with a flux of v electrons per 

second, the complete equation of dynamical balance becomes 
dn n exp (-a n 1 / 3 ) 



where the numerical values of 0!andT o are given by Equations IX-1 andIX-6, 
respectively. In the case of v = 0, Equation X-l can be solved in closed form 
as given by Equation IX-9. In the case of v =/ 0, numerical methods have to 
be applied. However, an analysis can be made with the following asymptotic 
method. 

Expanding the first term of the right side of Equation X-l about a point n G - 
corresponding to the initial value of the defect density, Equation X-l becomes 
approximately ( X - 1 /t 0 ), and 

■jjp = - [(v-Xn G exp (-kn Q ^ 3 )]- X(l- |- n Q V 3 ) exp ( - kn Q V 3 )/i 

= A-B/j,, X-2 

where \i = n - n Q . We therefore obtain 



When B >0, the time development is a slow process, but when B>0 the time 
development is rapid. 
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We have the following table, where the regions are shown in Figure X-l 
for case 2. 


Region 

A 

B 

Characteristic 

I 

+ 

+ 

Slow increase in damage 

II 

+ 

- 

Fast increase in damage 

in 

- 

+ 

Slow annealing 

IV 

- 

- 

Fast annealing 


The solid curve in Figure X-l represents n/r for various values of n where v 
is larger than a maximum of n/r . The four regions are represented by dif- 
ferent shaded areas. It is interesting to see that in region II, which represents 
a heavy initial damage, as time progresses, even at moderate irradiation rates 
the damage will increase further. 

In case 1, where v = max (n/r), regions IE and IV vanished. Therefore, no 
apparent annealing would be observed. On the other hand, if v is very small, 
regions m and IV will be expanded, that is, the annealing region will be 
dominant. 

The treatment here has been confined to the case where v does not change as a 
function of the time. In general, v is not a constant and the initial defect 
density has to be treated as a movable initial condition. The computation will 
be very complicated, but the principle involved is well illustrated by the 
present discussion. 
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XI. SOME EXPERIMENTAL RESULTS PERTAINING TO SPACE 
APPLICATIONS 

Some essential ingredients for simulating space conditions, from the point of 
view of radiation damage in solar cells, are: 

1. A spectrum of energetic rays and particles, mostly electrons and 
protons, with a moderate rate of about 10 11 to 10 13 particles/cm a -day. 

2. An ultrahigh vacuum. 

3. A system temperature from -100° to 50° C. 

In performing the experiments, for ordinary solar cells which have no signifi- 
cant annealing effects in the space environment, the radiation damage of solar 
cells depends only on the total radiation dose, but not an the dose rate. There- 
fore, the true function of the radiation v(t) is not important. As discussed in 
the section entitled "Dynamics of Annealing", in the case of the Li solar cell 
v(t) is of fundamental importance. The status of the present knowledge of some 
basic parameters required to predict the damage behavior for a given v(t) is 
such that even a semiquantitative estimation cannot be made. This implies the 
importance of obtaining information regarding the energy spectrum and dose 
rate of the radiation source. 

A high vacuum condition could be important and beneficial especially in the case 
of low energy irradiation. In a special experiment carried out in an ambient 
atmosphere, we observed that a prolonged radiation for several weeks on solar 
cells definitely leads to a discoloration of the silicon oxide antireflecting surface. 
This is caused by the surface bombardment by gas ions originating in the atmos- 
phere. 

The temperature effect of the radiation damage and annealing in the Li solar 
cell has not been studied sufficiently. Some initial observations imply that 
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the high damage introduction rate occurs at a temperature of -50° C or below. 

At temperatures higher than 70°C, a simultaneous annealing is so effective 
that no appreciable damage is observed. 

Some special experiments already carried out or now underway will now be 
discussed: 

1. Dose rate dependent experiment. In this experiment, van de Graaff 
irradiations of 1 Mev electrons of three flux rates are used, 1 x 10 14 , 1 x 10 16 , 
and 2.4 x 10 17 e/cm -day. The total accumulated flux is 3 x 10 14 e/ctoi a . 
Degradation of the short circuit current of three groups of solar cells, each 
group consisting of three, is measured 1 hour after the cessation of the 
radiation (Figure XI-1). 

First, we observe an increase of the residual damage as the flux rate is in- 
creased. An extrapolation indicates a vanishing residual damage if the flux 
rate is less than 1 x 10 13 e/cm a -day. This is generally satisfied in the actual 
space environment. 

Second, the data presented were for the residual damage measured 1 hour 
after the cessation of the radiation. A further annealing to complete re- 
covery is obtained if two additional hours are allowed to lapse in the 
1 XlO e/cm -day radiation case, while it requires 5 days for the 
2.4 xl0 17 e/cm 3 -day case. Therefore, a slow radiation rate results in less 
damage; also, slow radiation requires less time for recovery even when 
"time" is a sum of radiation time and annealing time. 
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The data of Figure XI-1 are obtained from solar cells fabricated 3 months ago. 
The experiment will be repeated with more recent solar cells in hopes of ob- 
taining better statistics from an improved type of solar cell. 

2. Repeated and interrupted radiation. In this experiment, a group of 
solar cells is irradiated with a flux rate of 1 x 10 16 e/cm a -hr to the desired 
level of accumulated flux (solid points in Figures XI-2, XI-3, and XI-4). The 
second group is irradiated at the same rate to 3 X 10 14 e/cm 3 * * * * (circular points 
in Figures XI-2, XI-3, and XI-4) and is repeated each time after about 10 days. 
The results show that the degradation in the case of repeated radiation and inter- 
ruption, for an equal accumulated flux, is smaller than in the case of a contin- 
uous radiation. In the same figure, a few triangular data points are given. 

This data represent the Type II Li cell which was described in the section en- 
titled "Different Types of Li Solar Cells". At high accumulated flux levels, 

the Typr II cell definitely shows a lower level of degradation. But in the case 
where the defect introduction rate is not excessive such a type of solar cells, 
which represents a somewhat complex construction, is probably not necessary. 

Although Figures XI-2, XI-3, and XI-4 give the statistical distribution of the 
experimental data, a simple arithmetic average, which gives a somewhat 
clearer picture, is given in Figures XI-5, XI-6, and XI-7. 

3. Some other experiments in prospect are: (a) Use a Sy 90 - Y 90 isotope 

source for /S-rays and irradiate at a level of 1 x 10 ia e/cm 8 -day in vacuum con- 

ditions at different temperatures; (b) determine the stability of the Li cell at 

different temperatures with different damaged states; and (c) determine the 

effect of photon injection on the rate of defect introduction and annealing. 
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XII. OUTLOOK 


Some background material and our present knowledge have been discussed 
so far. We now face the inevitable question of the comparison between the 
merits of the Li cell and those of the current n-on-p cell. In order to discuss 
this rather premature problem, some points to be kept in mind are enumer- 
ated as follows. 

1. Our knowledge of the structure of the defects as well as the 
mechanism of annealing is generally confined to some phenomenological 
observations. Therefore, we can predict neither the limit of the optimum 
efficiency in the initial performance nor the extent of the recovery which can 
be expected. This gap will be filled by an extensive program of basic studies 
centered on the following topics: Absorption band structure, local mode 
calculation, electron spin resonance, and Hall measurement of interstitial or 
displaced Li paired with radiation- induced defects in the same Si. Some of 
these programs are in progress and others are in the active planning stage. 

2. More quantitative measurements on solar cells will be made, with 
special emphasis on stability versus precipitation with regard to Li or pair- 
ing of Li with other defects. The flux dependence of the stability, especially 
at a very high flux level, will be a most important problem. 

3. Because of the observation of the dependence of the introduction rate 
as well as of the annealing on the defect concentration, analysis of the radiation 
damage requires a very specific description of the damage introduction rate 
and annealing condition. This includes the level of photon injection which 
produces a photovoltaic biased voltage; this in turn will influence the diffusion 
and the distribution of Li. 

4. Radiation damage due to different types of charged particle is in- 
teresting especially in the case of high-rate irradiation, as in the encounter 
with a solar flare, or in the case of neutron damage, where the defect density 
in some microscopic region could exceed that of Li. 
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In view of the above problems, it seems a correct conclusion that the im- 
portant question of actual application of the Li solar cell in space should be 
held in abeyance for a period of 6 months to 1 year until results are received 
from the extensive program which is now well in progress. 
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